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Abstract: We study the supersymmetric particle spectra and LHC collider observables 
for the large-volume string models with a fundamental scale of 10 11 GeV that arise in 
moduli-fixed string compactifications with branes and fluxes. The presence of magnetic 
fluxes on the brane world volume, required for chirality, perturb the soft terms away from 
those previously computed in the dilute- flux limit. We use the difference in high-scale gauge 
couplings to estimate the magnitude of this perturbation and study the potential effects 
of the magnetic fluxes by generating many random spectra with the soft terms perturbed 
around the dilute flux limit. Even with a 40% variation in the high-scale soft terms the low- 
energy spectra take a clear and predictive form. The resulting spectra are broadly similar to 
those arising on the SPSla slope, but more degenerate. In their minimal version the models 
predict the ratios of gaugino masses to be Mi : M2 : M3 = (1.5 — 2) : 2 : 6, different to both 
mSUGRA and mirage mediation. Among the scalars, the squarks tend to be lighter and 
the sleptons heavier than for comparable mSUGRA models. We generate 10fb -1 of sample 
LHC data for the random spectra in order to study the range of collider phenomenology that 
can occur. We perform a detailed mass reconstruction on one example large-volume string 
model spectrum. lOOfb -1 of integrated luminosity is sufficient to discriminate the model 
from mSUGRA and aspects of the sparticle spectrum can be accurately reconstructed. 
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1. Introduction 



The imminent advent of the Large Hadron Collider (LHC) is an excellent motivation to 
develop techniques to relate high energy string compactifications to observable low-energy 
collider physics. The LHC will be an unprecedented probe of the terascale and of the physics 
that stabilises the electroweak hierarchy. If supersymmetry is discovered at the LHC, it 
will be necessary to connect the collider observables and the spectrum of superparticles to 
a more fundamental theory such as string theory. 

Supersymmetric phenomenology is the study not of supersymmetry but of supersym- 
metry breaking: the undetermined parameters of the minimal supersymmetric standard 
model (MSSM) are those associated with the soft supersymmetry breaking terms, such 
as gaugino masses, scalar masses or trilinear A-terms. There are over one hundred such 
parameters in a general phenomenological parameterisation of the MSSM. The many differ- 
ent possibilities for supersymmetric phenomenology are determined by the many different 
possibilities for the soft terms. However, high scale constructions such as string compact- 
ifications contain far fewer independent parameters and so could be expected to lead to 
distinctive patterns of soft terms. 
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Broken supersymmetry is a property of the vacuum of the theory. The study of super- 
symmetry breaking vacua in string theory therefore requires control over both the moduli 
potential and the quantum corrections that enter into it. The quantum corrections play 
a crucial role for the large- volume models [1, 2] that we study here. The first construc- 
tions of IIB string models with all moduli stabilised [3] involved unbroken supersymmetry, 
which had to be broken by hand by the addition of D3-branes, with a gravitino mass 
that could only be lowered below the Planck scale by fine-tuning. It was subsequently 
realised that with the inclusion of quantum corrections [4] the moduli potential can ad- 
mit non-supersymmetric minima [5] and even naturally generate a hierarchically small 
supersymmetry-breaking scale [1,2]. This is the large volume scenario that we will explore 
here. 1 

Nonetheless, it is still a long way from the mere existence of low-scale supersymmetry 
breaking to actual phenomenology, and turning any given string compactification into LHC 
collider observables requires many steps. We will assume gravity-mediated supersymme- 
try breaking as the most directly motivated scenario from string compactifications; the 
modifications for other proposals are straightforward. 

1. The first requirement is that the compactification moduli be stabilised within a con- 
trolled approximation. This is necessary to ensure that the compactification is a good 
vacuum solution of string theory. 

2. The second requirement is that supersymmetry is broken in the vacuum, so that 
soft terms can be generated. We also require the supersymmetry breaking to be 
hierarchically small, in order that the soft terms appear at the TeV scale rather than 
near the Planck scale. 

3. The third requirement is a visible sector containing the Standard Model gauge group, 
with an understanding of the non-renormalisable couplings between the visible sector 
and the hidden sector moduli that break supersymmetry. 

4. The fourth task is to combine these couplings with the hidden sector supersymmetry 
breaking to compute the visible sector soft supersymmetry breaking terms at the 
high-energy compactification scale. Here we also would like an understanding of why 
the soft terms generated do not lead to large CP violation or flavour-changing neutral 
currents. 

5. The fifth task is to run these soft terms down to the TeV-scale, in order to compute 
the physical sparticle spectrum. 

6. The sixth task is to put such a spectrum through event generators such as [7, 8] and 
detector simulators such as [9] in order to generate collider observables. 

7. Finally, we also want an estimate of the uncertainties arising from the above six steps. 

lr The possibility of generating low-scale supersymmetry breaking from the moduli-stabilising potential 
has also been studied in corners of string theory different from that of IIB with fluxes [6] . 
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The aim of this paper is to carry out this program almost in its entirety for a specific 
and well-motivated class of string compactifications, the large volume models developed in 
ref. [1]. These models arise within IIB flux compactifications with moduli stabilisation and 
are characterised by broken supersymmetry with an exponentially large compactification 
volume. This allows the natural generation of hierarchies between mass scales, an extremely 
desirable feature. The large volume V lowers both the string scale m s and the gravitino 
mass 7713/2 with respect to the Planck scale Mp, 

M P M P 

mS ~ W m3 /2~— • ( L1 ) 



V refers to a dimensionless quantity: the volume measured in powers of the string-length 
l s . The dimensionful volume of the compactification manifold is VI®. From (1.1), an 



intermediate string scale m s ~ 10 11 GeV gives rise to TeV-scale supersymmetry breaking. 
The phenomenological implications of these models have been studied in [1,2,10-12] where 
requirements 1-5 above have been addressed. In this article we will first review the most 
relevant results from those references and complete requirements 6 and 7 above in order to 
make direct contact with potential observables at the LHC. 

The organisation of this paper is as follows. In section || we review the large-volume 
models and the moduli stabilisation that generates the exponentially large volume. We also 
describe the computation of soft terms and explain why leading order flavour universality is 
assured, summarising the results of refs. [1,2,12,13]. We also describe how magnetic fluxes 
perturb the soft terms away from those computed in the dilute-flux limit. We estimate 
the magnitude of the flux perturbation from the non-universality of the high-scale gauge 
couplings. Such corrections will generate an uncertainty in the high-scale soft terms that 
will translate into an uncertainty in the low-energy spectra and observables. In section 



we examine the low-energy spectra arising from the soft terms considered in section 2.3 
by generating random soft terms perturbed about the dilute-flux limit. We describe the 
generic properties of the resulting low-energy spectra and compare with those arising in 
mSUGRA or mirage mediation models. In section || we study collider observables for these 
spectra. We use counting observables to scan the properties of the randomly generated 
spectra, and show that for a sample spectrum sparticle masses can be reconstructed. In 
section [| we present our conclusions. 

We have made an effort to make this article self-contained. A phenomenologically 
minded reader may wish to skip the formal details of section 2 and start in section 2.3. 



2. Large- Volume Models 

Large-volume models represent a class of string compactifications, with all moduli sta- 
bilised, in which quantum corrections to the scalar potential naturally lead to a exponen- 
tially large volume. They were first found in references [1,2]. They are robust against 
additional quantum corrections, such as those of [14,15]. A recent detailed study of this 

2 The omission is in the lack of an explicit magnetised brane configuration realising the MSSM gauge 
group; this we simply assume can be achieved. 
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robustness is [16]. They have been applied to obtain string theory inflation [17], natural 
QCD axions [18], the scale for neutrino masses [19] and to low energy supersymmetric phe- 
nomenology [10-13] in which they provide a natural hierarchy for supersymmetry breaking 
with the large volume leading to an intermediate string scale. A comprehensive review 
is [20]. We start by reviewing their construction and properties. 

2.1 Construction 

These models arise as a rather generic limit of flux compactifications of IIB string theory 
in the presence of D3 and D7 branes. In N = 1 supersymmetric IIB compactifications 
the Kahler potential and superpotential for the moduli = S,U a ,Ti take the standard 
form [3,4,21,22], 

= -2 In (v + -^— ) J -hi(i JnAOj - \n(S + S), (2.1) 
W(&) = J G 3 A n + ^ Ae _aiTi , (2.2) 

i 

where the dependence on the complex structure moduli U is encoded in the Calabi-Yau 
(3, 0) form $7. G3 corresponds to the three-form fluxes and is linear in the dilaton S. 
We have included the leading a' correction to the Kahler potential, which depends on 
£ = — (^(3)x(M) / (2tt) 3 with x(M) the Euler number of the Calabi-Yau manifold M. Large- 
volume models require M to have at least two Kahler moduli T,, one of which is a blow-up 
mode, as well as a negative Euler number, i.e. x(M) < 0. The simplest model is that 
of Pp 1169 j, which we use as our working example. For this the volume can be written 
as [1,23] ' 

v =wS' 2 - T - 2 )- (2 - 3) 

Tj, = Re(T{>) and t s = Re(T s ) denote big and small cycles. The geometry is analogous to 
that of a Swiss cheese: the cycle TJ, controls the volume ('the size of the cheese') and T s 
controls a blow-up cycle ('the size of the hole'). 
The J\f = 1 moduli scalar potential is 

V = e k (k^DiWD-jW -3\W\ 2 ^j , (2.4) 

where D^W = diW + (diK)W. Dropping terms sub-leading in V, this potential becomes 
^ K^D^WD^W \(a s A s ) 2 ^T s e- 2a ° r ° fiW a s A s T s e- a ^ v£\W \ 2 

$=S,U 9s V 



in the limit V>1, where A, /x, v are numerical constants. The first terms of ( |2.5[ ) stabilise 
the dilaton and complex structure moduli at DgW = DjjW = 0. The remaining terms 
stabilise the Kahler moduli. The non-perturbative terms in r s balance against the pertur- 
bative corrections in the volume, and it can be shown that at the minimum of the scalar 
potential [1] 

c 

V ~ Woe as ( t s ~ InV, 
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where Wo is the value of the flux superpotential at the minimum of S and U fields, and 
c ~ £ 2 / 3 is a numerical constant. The resulting volume is exponentially large with a 
small blow-up cycle. This geometry is shown in figure p]. The minimum avoids two 




Figure 1: The physical picture: Standard Model matter is supported on a small blow-up cycle 
located within the bulk of a very large Calabi-Yau. The volume of the Calabi-Yau sets the gravitino 
mass and is responsible for the weak/Planck hierarchy. 



problems of the KKLT scenario. First, as its consistency does not require a fine-tuning of 
Wo, its existence is more generic. Secondly, the stabilisation of the T moduli can never 
destabilise the U moduli, as near the minimum the relevant terms in (2.5) come with 
different powers of the volume. A phenomenological advantage is that this minimum also 
breaks supersymmetry and generates an exponentially large volume. 3 

This large volume allows the generation of hierarchies by lowering both the string scale 
and gravitino mass as in (IT). A volume V ~ 10 15 is required to explain the weak/Planck 



3 The minimum of the potential (2.5) is at negative vacuum energy |Vo| ~ m^j^Mp with supersymmetry 
broken. Just as in KKLT a lifting term is desirable to obtain de Sitter or Minkowski space. There are 
several possible sources for this lifting term [24] (see [25] for a recent detailed analysis for KKLT and large 
volume models). Since the anti de Sitter minimum is already non supersymmetric, in contrast to KKLT 
models, the contribution of this lifting term to the soft breaking terms is suppressed and does not play an 
important role in the rest of this paper. 
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hierarchy and give TeV-scale super symmetry breaking. We assume such a volume through- 
out this paper. Further attractive features are that the volume V ~ 10 15 also gives an axion 
decay constant f a ~ 10 11 GeV in the allowed window [18] and the required neutrino sup- 
pression scale of 10 14 GeV [19]. 

A fully realistic model must have a Standard Model sector, which requires an ap- 
propriate configuration of O-planes and magnetised D7-branes. Chirality arises from the 
topological intersection numbers of differently magnetised branes. 4 To avoid having gauge 
groups that are too weakly coupled, we must wrap the Standard Model D7 branes on the 
small cycle corresponding to t s . As the branes should wrap a blow-up cycle, the brane 
configuration will represent a local construction of the Standard Model. We assume such a 
configuration can be found, but do not attempt to realise the brane configuration explicitly. 
The techniques involved in explicitly constructing such a configuration will be similar to 
those used in models of branes at singularities. In this context there has been recent effort 
at constructing MSSM-like gauge groups - see [26-28] for progress. 

As m s 3> ?7J3/2) the infrared physics of these models is that of the MSSM. The com- 
putation of the soft terms is then central to the study of the low energy phenomenology. 

2.2 Gauge Couplings 

The gauge kinetic functions f a (&) axe principally determined by the cycles wrapped by the 
D7 branes. If Tj is the Kahler modulus corresponding to a particular 4-cycle, reduction of 
the DBI action for an unmagnetised brane wrapped on that cycle gives 5 

U = (2-6) 

To generate chirality we require the branes to be magnetised. The magnetic fluxes alter 
dU) to 

fi = ^ + h l (F)S, (2.7) 

where hi is a topological function of the fluxes present on the brane. This can be understood 
microscopically. The gauge coupling 1/g 2 = Re/ is given by an integral over the cycle £ 
wrapped by the brane, 

.1= f dV~Vs + (2™')^. (2.8) 

Thus in the presence of flux, the gauge coupling depends on the magnetic flux as well as the 
cycle volume. 6 The factors hi(F) have been explicitly computed for toroidal orientifolds 
[29,30]. 

As a minimal scenario we assume that the Standard Model branes live on a single 
blow-up cycle and so the gauge kinetic functions depend upon only one Kahler modulus 



4 Magnetic flux on D7 branes is equivalent to dissolved D5 branes. The chirality arises from the point 
like intersection of the dissolved D5 branes. 

5 This holds with the phenomenology conventions Tr(T a T b ) = ±<5 ab ; with conventions Tr(T a T 6 ) = 5 ab 
the gauge kinetic function is / = 

6 In computing explicit expressions for the factors hi(F), it is easiest to use the Chern-Simons action to 
get the correction to Im(/), and then use holomorphy for the correction to Re(T). 
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T s . As the cycle size T is increased, the magnetic fluxes are diluted and their contribution 
to the gauge couplings goes away. It is useful to consider this dilute flux limit in addition 
to the physical case, and to see the latter as a perturbation on the former. The Standard 
Model gauge kinetic functions will be 



1. Dilute flux limit 



T 

fsu(3) = ^ 
T 

fsU(2) = ^, 

T 

fu(i) Y = k y^- ( 2 - 



2. Physical case 

fsu(3) = ^ + h su(3)(F)S, (2.10) 

fsu(2) = ^ + h su{2) (F)S, (2.11) 

fum Y =ky(^+h u{1) (F)sY (2.12) 

The factor ky accounts for the uncertainty in normalising the U(l) gauge field: the 
U(l) = U(n)/SU(n) that appears in intersecting brane models in general has a different 
normalisation to that of U(l)y. The 'canonical' value ky = 5/3 only holds for SU(5) 
grand unified models; in general, ky is model dependent. Typically in D-brane models 
hypercharge is an anomaly free linear combination of different [7(l)'s coming from different 
U (n) factors, the particular linear combination determining the value of the normalisation 
factor ky. The functions h(F) depend on the microscopic configuration of branes and 
fluxes, and we also regard these as unknown. 

2.3 Soft Breaking Terms 

A four dimensional N = 1 supergravity Lagrangian is specified at two space-time deriva- 
tives by the Kahler potential K, superpotential W and gauge kinetic function f a . The 
computation of soft terms starts by expanding these as a power series in the matter fields, 

W = W{§) + v($)HiH 2 + ^Y a ^)C a C^C^ + ..., (2.13) 

K = K(Q, <§) + K a -^, §)C a CP + [Z($, $)#i# 2 + h.c] + ..., (2.14) 

fa = /«(*)■ (2.15) 

Here $ denotes a generic modulus field and C a a generic matter field. In the MSSM the 
\i and Z terms apply only to the Higgs fields and so for these we have written Hi and 
H 2 explicitly. Gravity-mediated supersymmetry breaking is then quantified through the 
moduli F-terms, given by 

pm = e K/2K mn D n \fr. (2.16) 

The relationships of the expansions ( p.!3|) to ( 2.15| ) to the soft supersymmetry breaking 
terms are given in full detail in ref. [31]. 
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2.3.1 Gaugino Masses 

The canonically normalised gaugino masses are given by 

i F m d m f a 



M n . 



2 Re/ a 



(2.17) 



The gaugino masses that follow from (2/7) are 
1. Dilute flux limit 



2. Physical case 



Mi = M 2 



Mi 
M 2 
M 3 



F s 

— = M. 

2r, 



(2.18) 



2(r fl + 47r/ii(F)Re(5))' 
F s 

2(r s + 4 7 r/ l2 (F)Re(5))' 
2(t s + 4^Th 3 (F)Re(S)) , 



(2.19) 



where we write Mj for Msjj^i), and similarly for /i a and / a . 

In the limit of large cycle volume, the flux becomes dilute and the gaugino masses 
become universal at the compactification scale where the soft parameters are computed. 
In the physical case, the gaugino masses are non-universal due to the flux contribution. 
However for non-abelian gauge groups - i.e. for the wino and gluino - it follows from ( 2.1C] ) 
and (2.11) that the fractional non-universality of gaugino masses is identical to that of the 
gauge couplings: 

M 3 



M 2 



9l 



(2.20) 



Due to the factor ky in the U(l) gauge couplings, this relation does not hold for M\. Here 
we have 

M 3 



Mi 



si 



(2.21) 



where ky is the unknown normalisation factor. 
2.3.2 Scalar Soft Terms 

For the case of diagonal matter field metrics, K a g = K a 5 a p (no summation over a), the 
scalar masses, A-terms and B-term are given by [31] 



m 



1 - (mL +V )- F m F n drnd n log K a 



A 



Qtj3~f 



K m + d m log Y af3y - d m \og(K a KpK^ 



(2.22) 
(2.23) 



By. = (K Hl K H , 2 



2 < e 



k ' 2 li ( F 



K m + d m log 11 -d m \og{K H iK H2 ) 
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+ (2m\ /2 + V )Z- m 3/2 F m d m Z + m 3/2 F m [d m Z - Zd m \og(K Hl K H2 



pmpr, 



d m d n Z - {drnZ)d n log(K Hl K H2 ) }■ (2.24) 



As mentioned above, the tree-level vacuum energy Vo (in Planck units) is much smaller 
than (Vo/Mp < rriy 2 /Mp <C m 3/ 2 ) an d thus the effects of uplifting can be neglected. 
To compute the scalar masses and A-terms, the key piece of information required is the 
modular dependence of the kinetic terms K a . In the dilute flux limit this can be derived by 
relating the modular scaling of K a to that of the physical Yukawa couplings Y a ^ through 
the relation 



e 



K/2y 



Y^ = -~ (2-25) 



{K a K p K. 



i) 2 



and the fact that the T moduli cannot appear perturbatively in Y a ^ due to the combination 
of holomorphy and the Peccei-Quinn (PQ) shift symmetry. Ref. [13] used this reasoning 
to derive 

T 1/3 

K a ~ ^kaW, (2.26) 

where 4> refers to the complex structure moduli. In the dilute flux limit this leads to scalar 
masses and A-terms given by 



m a 

A*fh = — ^ — = —M. (2.27) 
4M 

B = - — . (2.28) 

The superpotential /i-term can be shown to vanish due to scaling arguments [13]: in 
the limit that V — * do, any non-zero superpotential /i-term would generate a mass for 
the Higgs field above the string scale. A non fine-tuned (i.e. weak-scale) /i term can be 
generated through the Giudice-Masiero mechanism [32], and as such is guaranteed to be 
of a similar order to the other soft terms. In fact it turns out to be difficult to satisfy the 



1 F s 


M 


7^2^ 


"7! 


F s 






—M. 






AM 




3~' 





B-term constraint (2.28) when running to low energies. However, the Higgs sector is by far 
the least understood sector of the Standard Model. As such it may be non-minimal, or the 
B-term may receive an extra contribution through the vacuum expectation value of a gauge- 
invariant scalar, a{N)HiH2- In this paper we therefore follow normal phenomenological 
practice by trading B for tan f3 and treating tan (3 as a free parameter. 

However, in addition to modifying the gauge couplings, the magnetic fluxes will also 
modify the kinetic terms and their dependence on t s . This affects soft breaking terms 
except in the limit of large cycle volume and dilute fluxes, where flux dependence should 
disappear. While there do not exist explicit formulae for Calabi-Yau compactifications, for 
compactifications on toroidal orientifolds such behaviour has been analysed using string 
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scattering techniques. A typical result is (see [29] and section 4.4 of [33]) 

K cc = {S + S)~ a IJ(T* + + Wr^^^^j^, (2.29) 

where Uj,S are the complex structure moduli and dilaton fields respectively, a, (3, 7 are 
constants and fa* = arctan (^T^J ^ s the relative angle between the branes, with p a flux 

quantum number. We can see that the T-dependence involves the angles fa. However, 
in the limit that T — * 00, fa — > and the flux-dependence disappears. We also see that 
fields with different fa, corresponding to different gauge charges, experience a different In- 
dependence. While toroidal examples are not the case of direct interest, they are important 
because they do allow an explicit computation of flux effects and one can explicitly see how 
the effects of fluxes become less important at large cycle volume. 

The large-volume models rely on a Calabi-Yau geometry and there do not exist any 
direct formulae for the effects of magnetic fluxes on the matter metrics. Nonetheless, 
the fluxes will affect the soft terms and must be taken into account. To model the flux 
corrections, we shall use the simple ansatz 

~ _ fo + e^F)) 1 / 3 

K <* ~ ylfz • ( 2 - 3 °) 

e a is used to parametrise the flux effects, and in the dilute flux limit e a <C t s . While the 
form of fl2.30p is simpler than will actually occur, it satisfies the basic requirement that the 
flux contribution will vanish in the limit that the cycle size goes to infinity and the fluxes 
dilute away. Using ([2.30D we can then compute 



V3 2(r s +e a (F))' 



A*/3 7 = --y= (m a + mp + m 7 ) . (2.31) 



In the limit that e a — ► 0, this recovers the dilute- flux expressions ( 2.27| ). 
2.4 Flavour and CP issues 



The formulae fl2jp to (gTgg ) used ab ove in computing the soft terms are a restricted form 



applicable for flavour-diagonal soft masses. This requires some explanation, as 'generic' 
gravity-mediated models give flavour non-universal soft terms and corresponding problems 
with flavour-changing neutral currents. The problem with this 'generic' expectation is that 
it is essentially a naive argument using only effective field theory, which does not take 
into account the actual structures that arise in string compactifications, which violate the 
genericity assumptions. 

It was argued in [12] that for the large-volume models, and more generally for models 
arising from IIB flux compactifications, there is a clean understanding of why supersym- 
metry breaking should give universal soft terms, at least at leading order. We briefly 
summarise the argument. The ability to distinguish flavours and consider flavour mixing 
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comes from the structure of the Yukawa couplings. The Yukawa couplings are generated 
from the superpotential, and as such can only depend on the dilaton and complex structure 
moduli. The combination of holomorphy and the PQ shift symmetry Im(T) —* Im(T) + e 
implies that the Kahler moduli cannot make any perturbative appearance in the super- 
potential: since flavour is generated in the superpotential, the interactions of the T-fields 
are flavour-blind. The physical Yukawa couplings also depend on the Kahler metrics; the 
scaling arguments entering ( [2.26 ) however guarantee that different flavours have the same 
scaling with the T-moduli. 

However, in IIB flux models it is the T-fields that have non-vanishing F-terms and 



break supersymmetry. The Kahler potential (2.1) also has a block-diagonal structure 



(K 



K, 



ss 



\ 








K, 



uu 








+ o 



K. 



(2.32) 



TT 



Mixing between the T and S, [/-fields is volume-suppressed and tiny. Thus the fields that 
break supersymmetry - the Kahler moduli - and the fields that give flavour - the S and 
U fields - are decoupled, to leading order, and supersymmetry breaking generates flavour 
universal soft terms. 

Large CP-violating phases are likewise not a problem. From the structure of the soft 
terms (2.17), ( 2.22; ) and ( p. 23 ), we see that the gaugino and A-term phases are all inherited 
from the small modulus F-term F s . They are thus universal and do not generate large CP 
violating phases that would be in conflict with observations. 



3. Spectra 
3.1 Generation 

The low-energy mass spectrum is determined by evolving the soft terms from the high scale 
to the TeV scale. Ambiguities in the high-scale soft terms will translate into ambiguities 
in the physical spectrum and observables. In ref. [12], a basic phenomenological analysis 
was carried out using the soft terms in the dilute- flux limit, ( 2.18| ) and ( p. 27 ), which were 



run down to produce TeV mass spectra. However, as emphasised above, the contribution 
of magnetic fluxes automatically introduces a theoretical uncertainty in the high-scale soft 
parameters. We want to understand how this uncertainty manifests itself in the possible 
low energy spectra. 

We make the assumption that the cycles are sufficiently large that the effect of the 
fluxes is as a perturbation on the dilute-flux results. In general the fluxes are unknown, 
but in one instance they can be taken 'from data'. In the dilute-flux limit, from fl2.9| ) the 
non-abelian gauge couplings would be exactly universal at the high scale m s ~ 10 11 GeV. 
From (|2.10| ) and ( 2.11 ), it follows that the fluxes are directly responsible for the difference 



of SU (2) and SU(3) gauge couplings. The magnitude of the flux perturbation can then be 
estimated from this ratio. 
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Assuming an MSSM spectrum with no exotic matter, the ratio of high-scale gauge 
couplings is 

9i 



9l 



« 1.37. (3.1) 

10 n GeV 

We regard both g| and g\ as fluctuations around a common value. In determining the soft 
terms that we run down to generate low energy spectra, we use the following strategy. We 



first specify a high scale value for the gluino mass, M3. The relations (|2.20| ) and (3.1) then 
fix the high-scale wino mass, 



M 2 



=1011 GeV 1.37 



' M 3 



,ou GeV (3 ' 2) 



The gaugino masses are flux- induced perturbations about a base value Jp-. We use the 
central value of M2 and M3 to estimate this, 

-L f» M c = — -. (3.3) 

We generate the remaining soft masses as fluctuations about M c . For example, 

Mt = M c (l±ei), 

m a = -|(l±e a ). (3.4) 

Here m a stands for the soft breaking mass for scalars. The perturbation parameter a a 
differs for each type of field, but is assumed to be the same across generations. 7 

With the ansatz (2.30), the A-terms are given in terms of the scalar masses by ( 2.31| ). 



A a (3j = --^{m a + mp + m 7 ). (3.5) 



We first generate a set of high-scale soft terms according to the relations ( |3.4[ ). For 
each spectrum, the e a were randomly generated within a domain < e a < eo with constant 
probability density. We initially take eo = 0.2, but also investigate the choice eo = 0.4. In 
generating spectra, as stated above we do not impose the high-scale value for the B-term 



( 2.28 ) and instead treat tan/3 as a free parameter, which we allow to lie in the region 
5 < tan /3 < 40. We repeatedly generate many random spectra in this manner. We then 
remove any spectra which fail experimental constraints, although all points pass direct 
sparticle search limits due to the heavy SUSY breaking scale set. 

Using the program S0FTSUSY2 . [34] the 2-loop renormalisation group equations (RGEs) 
for the MSSM are solved numerically to obtain a particle spectrum at the weak scale. The 
values of the Standard Model input parameters used for our computations are mt = 171.4 
GeV [35], m 6 (m t ) M = 4.25 GeV, a s (M z )^ = 0.1187, a" 1 ^) 3 " = 127.918, M z = 



7 This is required from considerations of flavour physics. The theoretical justification for this is that the 
flux magnitudes are dual to brane intersection angles 8i, and fields of different flavour but with the same 
gauge charges see the same angles 9i. It would be useful to further examine this question within explicit 
models. 
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91.1187 GeV [36]. Every particle spectrum generated must be regarded as equally consis- 
tent with the large- volume scenario we study. In determining what counts as an acceptable 
spectrum, we impose experimental constraints on the magnitude of BR{b — * 57), the 
anomalous magnetic moment of the muon (g — 2)„ and the dark matter relic density fth 2 . 
We require spectra to generate values for these within 2a of the experimental results. 

The average measurement of the b — > 57 branching ratio was obtained from [37]: 
BR(b — ► 57) = (3.55 ±0.26) x 10 -4 . For an estimate of the theoretical uncertainty we used 
the result of [38] 0.30 x 10" 4 ; adding the two errors in quadrature, we obtain the la bound, 
BR(b 57) = (3.55 ± 0.40) x 10~ 4 . 

The anomalous muon magnetic moment = (g^ — 2)/2 represents the largest the- 
ory/experiment discrepancy in precision electroweak physics. The average experimental 
value of is 116592080(63) x 10~ n [39], with an error that is statistics dominated. The 
dominant uncertainties in the Standard Model computation of arise from hadronic light- 
by-light scattering and vacuum polarisation diagrams - for recent reviews see [40,41]. The 
evaluation of the vacuum polarisation diagrams is carried out using experimental data from 
e+e~ -> hadrons. The Standard Model result given in [40] is af e = 116591785(61) x 10~ u , 
giving a 3.2 a discrepancy 

5a^ = (287 ±91) x 10 -11 . 

It is important to note that 5a^ usually has the same sign as fi, so a supersymmetric 
explanation of the — 2 discrepancy prefers \i > 0. 

We also impose bounds on the Higgs mass obtained by the LEP2 collaborations [42]. 
The lower bound is 114.4 GeV at the 95% CL. The theoretical computation of the Higgs 
mass in supersymmetric scenarios is subject to an estimated error of 3 GeV, and so we 
impose > 111 GeV on the result obtained from S0FTSUSY. 

The 2a WMAP [43] constraint on the relic density of dark matter particles is 

0.085 < Vth 2 < 0.125. (3.6) 

Assuming a thermal relic abundance, we compute the neutralino contribution to dark 
matter using mi cr Omegas 1 . 3 [44]. The lower bound in ( |3.6[) is only applicable if we require 
that the dark matter is solely composed of neutralinos. Since there may be other dark 
matter constituents, such as axions or for the large volume models the volume modulus 8 , 
when considering spectra we only impose the upper bound on Qh 2 . 

In determining the collider phenomenology the most important feature of any SUSY 
spectrum is the overall scale, in particular the scale of the squarks and gluinos whose pair 
production initiates the majority of supersymmetric events at a hadron collider. This has 
a large but mostly trivial effect on the observables, primarily through an overall rescaling 
of the number and energy scale of SUSY events: the lighter the spectrum, the more events 
that are generated. In studying the different spectra produced by the large-volume models, 
and how the high-scale flux uncertainty translates into a low-scale spectrum uncertainty, 

8 The volume modulus is light and has the potential to overdose the universe. Its abundance must 
therefore be diluted. We refer interested readers to [45] for a detailed discussion on this issue, as well as 
other astrophysical and cosmological implications. 
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our interest is not so much in the overall scale of the spectrum as in its structure. We 
therefore take a fixed value M3 = 500 GeV at the high scale, which corresponds to a 
physical gluino mass rrig ~ 900 GeV. Assuming supersymmetry is discovered at the LHC, 
the overall production scale could be constrained by a variable such as M e j t [46] . In section 
4 we will consider the effects of varying the overall scale. 

3.2 Features 

The particle masses of 200 spectra passing the 2a experimental constraints and randomly 
generated with M3 = 500 GeV at m s = 10 11 GeV and 20% variation parameters e a , are 
plotted in figure I 9 . From fi gure H we can see the extent to which uncertainties at the 
high scale translate into uncertainties at the low scale. It should be noted that because the 
fluctuations are assumed to be the same for each generation, but to vary independently for 
each type of field, there are correlations among fields of different generations. 

All of these spectra satisfy the above observational constraints. The tan (5 parameter 
was allowed to be in the range 5 < tan (3 < 40, but all of the points that pass the constraints 
have tan (5 < 23. The generic features of the spectra are 

1. The first-generation squarks are heavy, with masses around 800GeV. The : m g 
mass ratio is well-predicted, with a range from 780:900 to 810:900. The lighter stop 
(particle 12 in the figure) is at around 600 GeV, while the sbottom (particle 17) has 
a 720 GeV mass. 

2. The sleptons are much lighter at around 300 GeV. The fluctuations in the mj : m g 
ratio are much larger than for the squark and inherit their magnitude from the high- 
scale fluctuations. The stau is the lightest slepton, but is comparable in mass to the 
e and ft. 

3. The chargino X\ nas m ^ ± = 304 ± 5 GeV and exhibit very little fluctuation. It is 
nearly degenerate with the second neutralino x<]> which tends to be mostly wino. 

4. The LSP tends to be mostly bino, with a mass that can fluctuate substantially be- 
tween 200GeV and 300GeV. If the LSP has a mass towards the top end of this range 
it can have a substantial wino component. 

5. The charged Higgs fields are intermediate between the squarks and sleptons, with 
masses around 500 GeV. 

The choice of 20% for the variation parameter eo is somewhat arbitrary. We have also 
generated spectra in which the variations were only constrained to be within 40% of the 
central value M c . These spectra are shown in figure ||] 10 . 

9 The dilute flux spectrum shown is chosen to give a similar sparticle production scale compared to 
the other spectra shown in the figure. Another useful comparison would be a dilute flux spectrum with 
M3 ~ 430 GeV, which corresponds to the base value M c used for generation of the large volume models. 
This spectrum is not expected to lie at the centre of the large volume model spectra because of the different 
gaugino masses defined at the string scale, which in turn implies different RG evolution. 
10 See footnote H for comments on the dilute flux spectrum. 
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Figure 2: Particle spectra with the gluino mass normalised to ~ 900 GeV and 20% high scale 
uncertainties in other parameters. For comparative purposes, an SPSla spectrum with the same 
value of the gluino mass is also shown using black dots. The hollow circles show a spectrum coming 
from soft terms in the dilute flux limit. Particles 1-6 are the left and right handed sleptons, particles 
7-18 are the squarks, 19 is the lightest supersymmetric particle (LSP) Xi> 20-22 are the remaining 
Xi, 23-26 the Higgs particles, 27 the gluino, 28-30 the sneutrinos, 31 and 32 are the charginos Xi a- 



The spectra with 40% variation in high-scale parameters exhibits the same basic struc- 
ture as that with 20% variation, except that the spread is greater. The variation in squark 
masses remains much less than that for weakly interacting sparticles. 

The features of the spectrum described above can be explained analytically. The 
simplest example is for the ratios of gaugino masses. The one- loop RGEs for the gauge 
couplings and gaugino masses are 



dg a 
dt 
dMg 
dt 



b a , 



16vr 2 
16vr 2 



It follows from ( |3"77| ) and ( |3lj| ) that at one-loop 



and so 



Ms 
M 2 



M z 



9j 
9l 



d 
dt 



M z 



Ma 
9l 
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Figure 3: Particle spectra with the gluino mass normalised to ~ 900 GeV and 40% high scale 
uncertainties in other parameters. For comparative purposes, an SPSla spectrum with the same 
value of the gluino mass is also shown using black dots. The hollow circles show a spectrum coming 
from soft terms in the dilute flux limit. Particles 1-6 are the left and right handed sleptons, particles 
7-18 are the squarks, 19 is the LSP \ii 20-22 are the remaining Xi> 23-26 the Higgs particles, 27 
the gluino, 28-30 the sneutrinos, 31 and 32 the charginos Xi a- 



where we have used the high-scale gaugino mass expressions ( |2. 19 ). For the bino mass, we 
have 

M 2 



Mi 



Ah. 



9l 



M z 



\M ig 2 J 



M x 



kyg( 



(3.11) 



M z 



The bino mass cannot then be taken directly from data as it depends on the factor ky 
which enters the U(l) normalisation. The fluctuations in the bino mass at the weak scale 
are inherited from its high-scale fluctuations and are large. The bino tends to be heavier 
than in mSUGRA models. This can be understood from the leading order universality of 
the gaugino masses at the intermediate scale, which holds up to the effects of magnetic 
fluxes and tends to compress the gaugino mass spectrum. 11 It may be possible to achieve 
an M\ : M% ratio comparable to that of mSUGRA, but as seen in figure || it would be non- 
generic and difficult to achieve while still treating the magnetic fluxes as perturbations. 

To explain the structure of scalar masses and the magnitude of their fluctuations, 
we can use the approximate analytic solution for the soft parameters that is obtained by 



11 A similar compression is observed in mirage mediation models, where the gaugino masses are exactly 
universal at the intermediate scale. 
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one- loop RGE running to the low scale [47]: 
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(3.12) 



Here Ki is determined by the RGE running of the gaugino mass Mi, and is given by 



K a (Q) 




dtg 2 a (t)\M a (t)\ 2 (a = 1,2,3) 



(3.13) 



where Qo is the high scale and Q is the scale at which the squark and slepton masses are 
evaluated. The A contributions come from the D-terms and are expected to be small, 
being proportional to m|. 

For the large- volume models we numerically obtain K% ~ 1.92 M^(Mx) with M 3 (Mx) = 
500 GeV, Mx = 10 11 GeV. For comparison, an mSUGRA point on the SPSla slope has 
K 3 w 3.84 Mi(M x ) with M 3 (M X ) = 500 GeV, M x ~ 10 16 GeV. 

The effect of ( (Up and ( gl^) is that as we run to low energies the gluino mass 
comes to give the dominant contribution to the squark masses. This is due to the large 
contribution of the K 3 factor in equations ( 3.12j ). For small rriq(Mx), mq{MsusY) is mainly 
determined by the gluino mass and is relatively insensitive to the initial value of rriq(Mx)- 
Here Msusy is the scale of the SUSY breaking spectrum, ~ 0(1) TeV. This explains 
why the fluctuations in rriq(MsusY) are small, since the gluino mass is fixed at tree-level 
by M 3 = 500 GeV. The slepton evolution is however driven only by the wino and bino, 
corresponding to the weakly coupled SU(2) and U(l) gauge groups. RGE effects from the 
winos and binos are less pronounced than those of the gluino, and thus mj(Mx) gives the 
dominant contribution to mj{MsusY)- The physical slepton masses inherit the high-scale 
fluctuations of mi{Mx) and this produces the large slepton mass fluctuations observed in 
figures ^ and ||[ 

For the spectra of figures § and |3[ the stau is lighter than the other sleptons, but the 
stau-slepton mass splitting is not as pronounced as mSUGRA models with the same tan (3. 
This is because there is less room for the stau to evolve when running from the intermediate 
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scale rather than the GUT scale, and because the larger value of Mi tends to compensate 
for the effect of the large Yukawa coupling in the stau RGEs. 

It is promising that even allowing for large 40% variations in the high-scale soft terms 
does not alter the overall structure of the low-energy spectrum. We now describe how the 
spectra of figure ^ and |3] can be distinguished from mSUGRA or mirage mediation models. 

3.3 Discrimination from mSUGRA models 

The spectra that appear in figures ^ and || have properties of 'typical' gravity-mediated 
scenarios and are broadly similar to those arising from mSUGRA models. However, there 
are several important differences which enable these to be discriminated. 

The most obvious is the gaugino mass ratio. In mSUGRA the ratios of all three gaugino 
masses are set by the gauge couplings, with low-scale values 



(M 3 : M 2 : Mi 



M? 



(<?l = 92 ■ 9l) , r -6:2:1. (3.14) 



M 



Here g\ = |<?y with the GUT normalisation of U(1)y- This ratio follows from gaugino 
mass universality at the GUT scale. From the spectrum of figure ^, we see that for the 
large-volume models we have 



(M3-.M2-.Mt) ={g 2 3 :g 2 2 :k Y g Y ) . ~6:2:(1.5-»2). (3.15) 



M 



M? 



( 3.15 ) thus gives both a distinct prediction for the m g : mass ratio together with the 
expectancy that the : ratio will differ from that in mSUGRA. We can understand 
the larger value of the bino mass relative to mSUGRA from the fact that gaugino masses 
are approximately universal at the intermediate scale, up to the effects of magnetic fluxes. 
To achieve the mSUGRA gaugino mass ratios in the large volume models requires a very 
large fluctuation in the high-scale bino mass away from the value expected in the dilute- 
flux approximation. This is non-generic and difficult to achieve while still regarding the 
magnetic fluxes as perturbations. An observation of the gaugino mass ratios (3.14) would 
thus disfavour the large-volume scenario. 

The use of gaugino mass ratios as observables has one distinct advantage, recently 
emphasised in ref. [48]. The matter content enters both the gaugino and gauge coupling 



RGEs in the same way (through b a ) and so the one-loop derivation of the ratios ( 3.14 ) and 



(3.15) is independent of any extra charged matter that may be present beyond the MSSM. 
String constructions generically contain extra vector-like matter that will affect the running 
of the gauge couplings. The point is that such matter will equally affect the running of the 



gaugino masses and gauge couplings, and so the low-energy ratios ( [3. 14 ) and ( 3.15| ) will 



be unaltered. In this respect these represent predictions that do not rely on the details 
of the matter spectrum. Another advantage of the gaugino mass ratios ( 3.15| ) is that the 



derivation of these soft terms shows that these ratios persist even beyond the dilute flux 
approximation, indeed up to the point where non-perturbative string corrections in e~ Ts 
would become important in evaluating the gauge kinetic functions. 

There is a caveat to be added concerning the above relationship of gaugino masses 



(3.15). Two factors enter into the determination of the ratio Mi : M 2 . The first is requiring 
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that the high-scale perturbations remain in the dilute-flux regime. As an extreme example, 
within the 20% and 40 % variations above it is never possible to get M\ : M2 as 1 : 10; 
such a result requires high-scale soft terms that violate the dilute-flux assumption. The 
second factor is that in determining the spectra that count as allowed in figures Q and || 
above, the upper bound on the WMAP relic density constraint (Qh 2 < 0.11) has been 
imposed. As the bino mass enters crucially into the evaluation of the relic density, this 
imposes a further selection cut on allowed values for the M\ : M2 ratio, in conjunction with 
the slepton masses. 

The first constraint (remaining with the dilute flux assumption) prefers an M\ : Mi 
ratio larger than found in mSUGRA, i.e. greater than 1 : 2. The second constraint 
strengthens this preference. For example in the 40% diagram, the spectra with the lowest 
values of M\ are also seen to have low values for the slepton masses (this allows e.g stau 
coannihilation to reduce the relic density to the WMAP allowed region). This requires a 
large downward flux-induced perturbation of the slepton masses in correlation with that 
of the bino mass, which is disfavoured. 

The spectra shown in figures [2] and |3], and the resulting gaugino mass ratio of 3.15 



are a convolution of the above two effects. Actually, it is expected that this discussion 
will need further modification. The reason is that the use of the thermal relic abundance 
computation relies on the assumption that cosmology is thermal from a temperature of T ~ 
O(10GeV) down. Such reheating temperatures are in general hard to achieve within string- 
derived theories due to the cosmological moduli problem. In the large-volume models, 
there exists a light gravitationally coupled O(MeV) volume modulus. Such a field has a 
lifescale longer than the universe and must be diluted for it not to spoil nucleosynthesis. 12 
Some dilution mechanism, such as thermal inflation, is therefore necessary to dilute the 
modulus. The presence of such a dilution mechanism alters the late-time cosmology. The 
proper computation of the relic abundance should start with the reheating temperature that 
applies after the dilution mechanism has been in operation. The dilution mechanism may 
itself lead to misalignment of the light modulus and thus require a reheating temperature 
much less than O(10)GeV. 

The upshot of this discussion is that to avoid the moduli problem the reheating tem- 
perature may end up being significantly less than assumed in the thermal computation 
of the relic abundance, and this will tend to modify the spectra and gaugino mass ratios 
above. A lower reheating temperature will imply the need for more efficient particle anni- 
hilation. In general this will tend to prefer greater degeneracies in the sparticle spectrum 
(e.g. to allow efficient coannihilation channels), but the detailed effects of a lower reheating 
temperature are beyond the scope of this paper. 

We note the spectra (|3.14|) and (|3.15|) could potentially be distinguished within the first 



year of LHC running; for example a di-lepton edge measurement m^o — m^o = 50 GeV to- 

A 2 Al 

gether with an estimated SUSY production scale of msusY ~ 1000 ±300 GeV is consistent 
with the second gaugino spectrum and not with the first. 

A more subtle distinction can be made through the relative values of the scalar masses. 



2 The cosmological properties of the large volume models are discussed at more length in [45]. 
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For example, the evolution of the first-generation squarks is largely driven by the gluino 
and the RGEs for these squarks exhibit focusing behaviour. If we neglect all but the strong 
dynamics, we have a solution of the one- loop RGEs [49] 

m ^ -V r 2) | r 2 ™|(M*) where T _ 9l(M X ) 
M|(^)-9 (1 T)+r Mj{M x ) ^ r " ' (3 - 16) 



As the renormalisation scale /x is lowered, an infra-red stable fixed point of r — > is 

approached, i.e. (^~M^J ~~ > §■ However, how close one gets to this infra-red fixed point 

limit for fi = Msusy depends upon the starting scale of evolution: for Mx ~ 10 16 GeV 
(i.e. for the mSUGRA case), r = 0.34, whereas for the large volume models, Mx ~ 10 11 
GeV, r = 0.45. Equation ( |3.16 ) shows that it is generally hard to obtain m g <C m g ; if 



we start with rriq ~ m g then this remains the case, and even if initially rriq <C m g then 
the squark will rapidly run up to obtain a mass rriq ~ m g . The point is that by starting 
the evolution at the intermediate scale, there is simply less time for the rriq : m g ratio to 
evolve compared to starting at the GUT scale. If the squark masses start below the gluino 
mass, as holds in our case, they have less time to evolve and so tend to be lighter than 
for an evolution commencing at Mqut'- the focusing behaviour of equations ( |3.16| ) is less 

efficient. The ratio (^^j for an intermediate scale model will then always be less than 

for a corresponding GUT-scale model with the same soft terms. 

This is manifest in our spectrum. We obtain a ratio rriq x : M g ~ 800 : 900, and which 
can be significantly smaller (down to : M g ~ 770 : 900 with the 40% variation). A 
similar choice of high-scale soft terms in mSUGRA gives : M g ~ 850 : 900, illustrating 
the greater running of the squark masses starting from the GUT scale. In the mSUGRA 
framework, lighter scalar masses at the GUT scale will reduce the : M g ratio at the low 
scale. However the lighter scalar masses also give rise to lighter sleptons. This can be seen 
by examining the SPSla spectrum in figures |2] and ||. While the squark masses are now 
relatively close to those occurring for the large volume models, at m gi : M g ~ 830 : 900, 
the sleptons are significantly lighter. It is thus not possible to fit the scalar spectrum of 
figures ||| and || with mSUGRA models unless the assumption of universal scalar mass at 
the GUT scale is relaxed: with mSUGRA either the squarks are heavier or the sleptons 
lighter than for the large volume models. 

This illustrates a general feature of the large-volume spectrum, which is that it is more 
compressed than those appearing in mSUGRA. This follows primarily from the fact that 
the soft term running starts at the intermediate scale rather than the GUT scale. The 
masses therefore have less time to separate compared to mSUGRA models, and thus the 
spectrum is more bunched. 

Clearly, the precise rriq : m g and mj : m g ratios are not quantities that can be rapidly 
measured at the LHC and would require years of high-luminosity running. However, while 
difficult to measure, these ratios are very interesting because they offer the possibility of 
estimating the scale at which the soft terms have been defined, and thus even the possibility 
of indirectly measuring the string scale. 



- 20 - 



3.4 Discrimination from Mirage Mediation Models 

A scenario that has recently attracted attention is mirage mediation [50-52] - see [53— 
56] for recent work. This corresponds to soft terms arising from supersymmetric KKLT 
stabilisation, with supersymmetry broken by an anti-D3 brane. The gravitino mass is 
determined by the fluxes and is naturally at the Planck scale. The hierarchy is generated 
by fine-tuning the fluxes to reduce the gravitino mass to a TeV. The soft terms arise from 
a combination of gravity and anomaly mediation. The soft terms are defined at the GUT 
scale and exhibit mirage unification at an intermediate scale, 



where a is the ratio of gravity to anomaly mediation and is usually taken to be 0(1) (see 
however the discussion in section 5 of [12]). In terms of this parameter the gaugino mass 
ratios are given by [48] 



The rriyy : mass ratio can be similar to that arising from the large- volume models, 
and thus it will not be possible to distinguish mirage mediation and large volume models 
from the : mass ratio. However the m g : ratio is substantially smaller for 
mirage mediation than for the large volume models, and thus here discrimination is possible. 
For example, assuming a ~ 1, a wino mass of 300 GeV would correspond in the mirage 
scenario to m g ~ 570 GeV and in the large volume scenario to m g ~ 900 GeV. While 
the gluino mass might be difficult to measure directly, in both models the gluino mass is 
correlated with the squark masses, which are easier to measure. Thus by measuring the 
m^o—m^o and mq L —m^o mass differences it will be possible to distinguish these two models. 
For example, suppose m^o —m^o was measured as 75 GeV, together with m^o > 200 GeV. 
Then a measurement of rriq L — m^o ~ 250 GeV would prefer mirage mediation models to 
large- volume models and a measurement m gh — m^o ~ 550 GeV would prefer large-volume 
models to those of mirage mediation. 

Compared to mirage mediation, the large volume models have a less bunched spectrum. 
This can be understood by the nature of the soft term universality that arises at the 
intermediate scale. In the large volume models, this is approximate and is broken by the 
magnetic fluxes on the brane world volumes. The effect of this flux-breaking is to raise 
the gluino mass in relation to the wino mass: as the gluino mass runs up, this broadens 
the low-energy spectrum. In mirage mediation models the gaugino masses exhibit exact 
universality at the intermediate scale, and so at low energies the m g : ratio remains 
smaller than for the large- volume models. 

It thus should be possible to use the gaugino masses to distinguish the soft terms 
produced by the large-volume models from those appearing in either mSUGRA or mirage 
mediation scenarios. Distinction from mSUGRA is possible through the rriyy : ratio 
but not through the m g : m^, ratio; distinction from mirage mediation is possible through 
the m g : ratio but not through the : ratio. In both cases further distinction 




(3.17) 



(M 3 : M 2 : Mi) 



(6 - 1.8a) : (2 + 0.2a) : (1 + 0.66a). 



(3.18) 
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should be possible through the spectrum of scalar masses. However it should be remarked 
that the mirage mediation scenario is also based on IIB string compactifications and the 
presence of magnetic fluxes may also affect the structure of soft terms in that scenario. 

4. Collider Observables 

In this section we investigate the collider phenomenology that follows from the spectra 
given in section ||. We analyse this through Monte Carlo simulation of 10fb _1 of LHC data 
for each spectrum. For one given model we simulate lOOfb -1 of LHC data and show how 
well the spectrum can be successfully reconstructed. 

4.1 Collider Model and Data Generation 

To generate events we used PYTHIA version 6.400 [8] linked to the PGS (Pretty Good 
Simulation) detector simulator [9]. We take the S0FTSUSY2.0 spectrum and link it to 
PYTHIA with the SUSY Les Houches Accord [57]. For each particle spectrum generated 
we have simulated 10fb -1 of data for pp collisions at 14 TeV. We also simulated the tt 
and WW/ZZ/WZ Standard Model backgrounds. We did not simulate the W/Z + jets 
background due to the large amounts of CPU time required. Data was originally generated 
in the LHC Olympics [58] format. This contains the particle energies and momenta for all 
particles (electrons, muons, hadronic taus, jets and photons) in an event, and b-tagging 
information for jets. The data analysis was performed using ROOT [59]. 

The SUSY production cross sections for a given spectrum can be computed using 
Prospino 2 . [60] at next-to-leading order in QCD. With a gluino scale m g ~ 900 GeV and 
the spectra of figure § and ||, the dominant production cross-sections are a gg ~ 0.53 — 0.54 
pb, (Tqg ~ 2.9 — 3.0 pb and Oqq ~ 1.4 — 1.5 pb. Squark-gluino production is therefore most 
significant, followed by squark-squark and then gluino- gluino production. 

Using PYTHIA and PGS, we simulated 10fb _1 of mock LHC data for each of the spec- 
tra that were generated in section ||. The basic cuts used were the PGS Level 2 triggers, 
which are summarised in the appendix. As all of these models have a similar SUSY pro- 
duction scale and cross-section, differences in the number of events passing the cuts must 
be attributed to the detailed structure of the spectrum rather than to the overall scale. 
To compare and contrast many different models, all with the same high-energy origin and 
with the same overall scale, we first use counting observables. These will probably not 
be very useful for discovery of physics beyond the Standard Model, being too sensitive 
to experimental systematics and physics unknowns such as parton densities and parton 
shower/matrix element approximations. Nevertheless, we may imagine a time in the fu- 
ture when these effects have been measured to some precision at the LHC after a beyond 
the Standard Model discovery, and counting observables might be used for model discrimi- 
nation. Their importance here is that they provide a simple and easily visualisable measure 
of the differences in the possible phenomenology across a wide range of models. Although 
not our emphasis in this study, the observables used are mostly chosen to be those in [61]. 13 

13 The results given for the large- volume models in [61] differ from those here as the soft terms used are 
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While the cuts employed here are somewhat different to those used in ref. [61], a rough 
comparison between their results and ours may be useful. We also investigate the effect of 
increasing the high scale parameter fluctuations as well as varying the overall scale of the 
soft terms. 



In section 4.2 we discuss the potential to reconstruct the spectra from kinematic ob- 



servables. We will find that this is much easier for some spectra than for others. 
4.2 Counting Observables 

A counting observable is simply the number of events in a sample which all satisfy a certain 
set of desired properties. Since the expected number of events N with properties V has a 
statistical uncertainty yN, the signal is taken to be observable only if it is well above a 
large statistical fluctuation of the background. Therefore for observability we require 

^ > 4, -f > 0.1, N s > 5, (4.1) 



N b N b 

where N s is the number of events in the signal data sample satisfying property V and Nj, 
is the number of events in the background sample satisfying the same property. Our esti- 
mation of Nf, could easily be wrong by factors of a few due to theoretical and experimental 
uncertainties, but we shall see that thankfully most observables will still be usable. 

The basic cuts for the counting observables are the L2 triggers used in the LHC 
Olympics version of PGS. On top of that we impose additional cuts: 

1. For all jets entering the counting, we require jet Pt > 100 GeV. 

2. For all isolated e, /i, we require Pt > 10 GeV. 

3. For all rs require Pt > 100 GeV. 

4. Missing transverse momentum j/t > 300 GeV. 

The choice of j/t cut can be motivated by the plot in figure |], produced for one of the 
randomly generated spectra. It can be seen that for pr > 340 GeV, the signal dominates 
over the background. This value does not depend significantly on the choice of the random 
spectrum. 

The simplest signatures to consider are inclusive counts of opposite sign di-lepton 
and tri-lepton events (including taus). These are also among the most promising for SUSY 
searches since the detector tagging efficiency for leptons (electrons and muons in particular) 
is quite high. Furthermore, there is relatively little Standard Model background for multi- 
lepton multi-jet high-yy events. 

The presence of many opposite sign di-leptons is also an indicator of whether the decay 
chain X2 — ► / ± / =F — ► l + l~Xi occurs frequently in the samples. This decay chain may be 



different: in [12] it was shown that a cancellation exists in the scalar soft terms that significantly reduced 
the scalar soft terms by a factor ln(Mp/m 3 / 2 ) compared to the original estimate of [11] that was used 
in [61]. The Kahler potential was calculated for Calabi-Yau manifolds to first order in a volume expansion 
only recently [13]. The cancellation occurs only for chiral fields. 
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Figure 4: j/t plot for background and a sample signal spectrum. The background is shown in red. 
The spike in the background j/t at 200 GeV is an artifact of the triggers used. 
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used [62] for constructing an l + l endpoint which can subsequently be used to constrain 
sparticle masses. The results for spectra with nig fixed at ~ 900 GeV are shown in Figure 

I 

Even though the overall mass scale of 
the spectrum is fixed and the SUSY cross- 
section is essentially unaltered between mod- 
els, there is still a large range for the num- 
ber of observed di-lepton and tri-lepton events. 
The dominant background comes from tt 
events, and in a large number of cases the 
number of such SUSY events lies above the 
acceptable number of background events. 

This variation can be understood in terms 

of the different possible spectra, The num- _ , T . ,, . . ,„„ ,. , 

Figure 5: Numbers of tri-leptons/OS di-leptons. 

ber of di-lepton (and hence tri-lepton) events The black dot corresponds to the observable limit 
depends crucially on the details of the X2 ~~ > according to ([O]). 
l^P — > I decay chain. In most cases 

the first 2 generations of left handed sleptons are heavier than ■ However the right handed 
sleptons may be lighter or heavier than yQ. If they are an heavier, then the 3-body decays 
of the X2 wm dominate. The X2 ~^ Xi decay will then often proceed through an off-shell 
slepton I*, giving many di-lepton events. However, if eR, flu are heavier than ^2 but f\ is 
lighter - which can occur as, for larger tan f3, is driven down by the RGEs and tends to 
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the tau tagging efficiency is relatively low, and nif — is relatively small (meaning that 
the taus tend to be soft), few taus will be picked up by the detector. It may also happen 
that all of en, fin, tr, v are lighter than yQ, with the dominant branching ratio of again 
into fs and vs. In this case we again observe only a few di-lepton or tri-lepton events. 

In ref. [61], the number of clean di- and tri- leptons was used as an observable, meaning 
no jet activity in the detector. From the point of view of our simulation, this corresponds 
to direct weak gaugino and/or slepton production. In 10 fb _1 , none of our model samples 
produced a single "clean" event, consistent with the masses of the sleptons/weak gauginos 
and the fact that this is a weakly interacting production channel. The usefulness of this 
observable is questionable, since there will always be some jet pollution in the detectors due 
to the QCD background. Thus some cut on hadronic activity must be given experimentally 
in order to define a jet veto, and the predicted backgrounds can be notoriously unreliable. 

Another observable considered in [61] that we will not use here concerns the number of 
events with no leptons, 1 or 2 6-jets, and at least six hard jets. The difficulty in using this 
observable comes in the estimation of the background. The processes in PYTHIA are 2^2 
rather than 2 —* many, and so the 'hard jet' background arising from PYTHIA comes from 
the parton shower rather than from direct hard jet production. This may underestimate 
the background by orders of magnitude. A correct estimate of this background would 
require the inclusion of 2 — > 4, 5, 6 processes in the Monte Carlo (e.g. with ALPGEN [63]). 

In figure |(| we plot the number of di- 
lepton multi-jet events with 1 or 2 6-jets 
against the number of di-lepton multi- 
jet events with 6-jets. There is a clear 
positive correlation in Fig. |(| indicating 
that the number of events with two lep- g 150 
tons and 6-jets is dependent more on the 
number of events with leptons passing 
the cuts than on the number of 6-jets. 



350 
300 
250 
200 



100 
50 



random LV, M3=500 



*w * * - 



Nonetheless the number of events with 100 200 300 400 500 eoo 700 soo 900 1000 

, . . i r p f-» 1 2 lep, b and >-2 jets 

o-jets can vary by a factor of 2 between 

different spectra. This is reasonable, be- Figure 6: Number of events with 2 leptons, 1 or 2 b- 
cause there are numerous competing de- jets and 2 jets/Number of events with 2 leptons, b- 
cays where the branching ratios are de- J ets and 2 j ets - The black dot in both case represents 
pendent on various independent input the observable limit according to Q. 
parameters. We give one example. If 

kinematically allowed, as here, the gluino has a significant branching ratio to tt\, which 
generates events with multi-jets and bs. The branching ratios of t\ — > ix? 2 an d bxt depend 
on tan/3 and At, since they affect the mixing of the stops. In figure ^ we also require two 
leptons. However the number of leptons observed depends on the mass differences between 
the light neutralinos and the sleptons. The point is that these parameters mentioned above 
are varied independently in our models, which explains why spectra with similar numbers 
of leptons can have different numbers of 6-jets. Other decay chains, for example g — > bb, 
can be analysed in a similar fashion. 
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Figure 7: Comparison of results with 20% and 40% variations: di- and tri-leptons counts and 
events with 6-jets. 



The number of 6-jets may potentially be used as a coarse discriminator between differ- 
ent high-scale constructions. Suppose the scalars are much heavier than the gluino, with 
a large Higgsino component to the LSP. In this situation, the gluino will primarily decay 
through a 3 body channel to bb~Xi 2 " as weu as other channels involving tops if kinemati- 
cally allowed - due to the large higgsino coupling to stops and sbottoms. These decays will 
all result in a large number of 6-jet events. An example of such a construction is the focus 
point region of mSUGRA. Since in our construction the light neutralinos are always gaug- 
ino dominated, the number of 6-jet events is expected to be smaller. We have verified that 
a focus point spectrum does indeed give many more 6-jet events than those arising from 
the large-volume models. Thus, once 6-jet tagging is understood, the number of 6-jets can 
potentially be used to distinguish different high-scale constructions with a similar overall 
production scale. 

Since the 20% random variation we 
chose for all parameters is an arbitrary 
choice, it is useful to compare the results 
with a sample of spectra with 40% varia- 
tions. The results (again with M3 = 500 
GeV at the high scale) are shown in fig- 
ures f?]-||. In figure [7| we see the number 
of di- and tri-lepton events, as well as the 
number of di-lepton events with or with- 
out 6-jets. The general structure of the 
observables is similar to that seen above 
with the 20% variations. However the 
40% variation spectra, although with the 
same overall mass scale (as defined by the 
gluino mass) can have significantly more 
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Figure 8: Number of events with leptons, 1 or 2 
6-jets and 2 jets vs. number of events with 2 leptons, 
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Figure 9: Comparison of event numbers with 20% variations and varying mass scale, and results 
for the SPSla slope: (a) numbers of tri-leptons/OS di-leptons, (b) number of events with 2 leptons, 
1 or 2 6-jets and 2 jets/number of events with 2 leptons, 6-jets and 2 jets. 

di-lepton events than the 20% variation spectra. This is due to the fact that in the 40% 
case, the LSP is more likely to acquire a significant wino component, in which case the 
large left-handed couplings result in a lot of lepton production Also it may happen that all 
the right handed sleptons and left handed sleptons are lighter than an d that the mass 
differences are all large, which again results in many observed di-lepton events. 

An interesting counting observable to consider is the number of events with one or two 
6-jets, requiring 2 hard jets in the event. In figure |8] we compare the number of di-lepton 
events without 6-jets against the number of lepton events with 6-jets. This observable 
is shown on the abscissa of figure |8|, and it is consistent throughout the entire sample of 
spectra. The observability limit defined by Eq. ( |4.1| ) is at (209,49) and is omitted from the 
figure. 

4.2.1 Varying the Sparticle Mass Scale 

We next investigate the effect of varying the overall mass scale which has so far been set 
by fixing the gluino mass. An arbitrary value for M3 is selected and 20% variations on the 
parameter F s /(t s ) allowed. The results are shown in figures |9]-|lC|, in conjunction with the 
results for the SPSla slope given by mo = OAmi/2, Ag = ~~ ^o^sgn// = +1, tan/3 = 10. In 
the plots, 50 SPSla points from M 12 = 250 GeV to M 12 = 887 GeV in steps of 13 GeV are 
shown. 

The main conclusion from the study of counting observables is that, even with the 
spectra restricted to the form of figures || and ^ the total number of observed SUSY events 
can vary widely. Even fixing the overall scale of the spectrum, the large-volume models 
still lead to a widely varying number of triggered events: models at the same scale and 
coming from the same high-scale theory can give quite different results as small changes in 
high scale parameters can lead to significant changes in observables. We do not make an 
explicit comparison of our results with the models discussed in [61] for this reason. If there 
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is so much variation within our models, there will be a larger variation with other models 
that would make them difficult to differentiate. We have also used different cuts, which 
makes a quantitative comparison impossible. Nonetheless, a rough look at the similar plots 
in both works does not indicate an easy way to separate these models from those in [61]. 



4.3 Potential for Reconstruction 

We now discuss the potential for recon- 
structing the spectra of figure ^. This 
will illustrate the above point: if there 
are very few di-lepton events, and no kine- 
matic endpoints, then direct reconstruc- 
tion would be very difficult if not impos- 
sible. The structure of the analysis given 
here follows standard accounts of recon- 
struction such as in ref. [62,64], for ex- 
ample. 

The ability to reconstruct supersym- 
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Figure 10: Comparison of event numbers with 20% 

variations and a varying mass scale. The results for 
metric particle masses depends significantly the spsia slope are included for comparison . 

on the spectrum and on the decay chains 

and their branching ratios. Since jet observables usually suffer from large combinatorial 
backgrounds, the cleanest measurements are those involving only leptons. Therefore the 
first step in reconstruction of a supersymmetric spectrum is a measurement of a di-lepton 
edge from the x[> — * Xil P chain. As explained in section O, whether we observe few or 



many di-lepton events from this decay chain is determined by the mass difference between 
Xi and x% and the slepton branching ratios. Figure |ll] shows the plot of the di-lepton 
invariant mass Mn = (pi x +Pi 2 )' J '(pi 1 +Pi 2 )/i f° r lOfb -1 of data for two spectra, one with 
many opposite sign, same flavour (OSSF) di-lepton events in the signal, and one with very 
few. There is no evidence for an edge in the di-lepton invariant mass in Fig. |Tl|b, If there 
are few di-lepton events, the spectrum is much harder to reconstruct since one has to resort 
to multi-jet observables. We therefore restrict to considering a spectrum which generates 
many OSSF di-lepton events, such that the di-lepton edge can be easily reconstructed. For 
this spectrum we simulated lOOfb -1 data and the plots given are based on this simulation. 
This represents one year of LHC running at design luminosity. The spectrum we attempt 
to reconstruct is as shown in Table [|. 

We start by selecting events that pass a set of cuts that we name selection A: 

1. tfr> 300 GeV 

2. Two opposite sign electrons or muons with Pt > 10 GeV. 

3. At least four jets with Pt 1(2)(3)C4) > 100(50)(50)(50) GeV. 

4. yf T > 0.2M e// , where M eff = P Tl + Pt 2 + Pt 3 + Pt a + Pt- 
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Figure 11: Opposite sign, same flavour di-lepton invariant mass for (a) a spectrum with many 
di-lcpton events and m^ 2 3> anc ^ O 3 ) a spectrum with few di-lepton events and m^o —m^o rs 21 
GeV. In both cases we require 4 hard jets and j/t > 200 GeV. Standard Model background is shown 
in red. 



m-± 

Al 


m xf 




m H o 




%± 










303 


480 


114 


532 


532 


538 










m~ g 


m d L 


mu L 


mu R 














909 


800 


792 


779 


583 


790 


725 


782 


787 




m e - L 


mfx 


m„ e 




m e R 


TOf 2 


Al 


A.2 


771^0 
A3 




348 


261 


338 


336 


270 


349 


233 


303 


460 


483 



Table 1: Mass spectrum of model picked for reconstruction. All masses are listed in GeV. The 
first two families are mass degenerate. 



We then plot the histogram of the di-lepton invariant mass, with a flavour subtraction 
of the e + fj>~ + e~/i + result to cancel processes with leptons arising from two independent 
decays. 

This plot is shown in figure |12[ It is well-known [62] that the xl ~> ^ T ~* Xi^F 
decay chain admits an endpoint at 



M max 



(m|. 



m~ )(Tn~ 



R Ai 



mi 



(4.2) 



The location of this endpoint can be found by fitting figure 12 with a triangular edge, 
smeared with a Gaussian (the width of which is also fitted) to simulate resolution effects 
of the experiment. MINUIT [66] and MINOS are used for this purpose, and to estimate 
the error on the measurement. We obtain M™ ax = 69.20 ± 0.15 GeV. 

Following Ref. [62], we next study the decay channel qi — > qx® ~^ l^l T q — * Xi^^Q to 
obtain a set of constraints on the squark masses as well as en, XijXE- Events are selected 
with the following properties (selection B): 
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Figure 12: Mu after cuts selection A. We expect an edge at ~ 69 GeV from the spectrum. The 
Standard Model background is shown in red. 



1. At least 4 hard jets, with P T ,i(2)(3)(4) > 100(50) (50) (50) GeV. 

2. Two opposite sign electrons or muons with Pt > 10 GeV. 

3. M eff = P Tl + P T , 2 + P n + P T4 +rf T > 400 GeV. 

4. rf T > max(300,0.2M e// ). 

The di-lepton 4- momentum is combined with each of the two hardest jets to obtain two 



different qll invariant masses. The lighter of these is plotted in figure 13. MJJ^ X can be 



written in terms of the maximum of several terms that contain sparticle masses and have a 



form similar to Eq. 4.2. In refs. [62,64], rather simple expressions were given for edges such 
as MJJ^ X . These were correct for the particular mass spectra examined in those papers, 
but are not true in the general case. We therefore use the general expressions given in 
Ref. [65] and refer the reader there for further details. An empirical fit of the form 

r M u ax / i \ 

f(M)=J " dz(a 1 (Ml[ l q ax -z)+a 2 (Mjr q ax -zf)e W ^-^ 

is used to reconstruct this endpoint with a = 25 GeV. With a variable a in the fit, MINOS 
did often not converge, indicating a possible degenerate x 2 minimum valley. With a fixed 



width, we were able to obtain a perfectly good fit, as Fig. |13| shows. The endpoint obtained 
through this fit is MJ^ X = 450 ± 5.5 ± 2.4 GeV. Where we quote two uncertainties, the 
first is a statistical one from the fitting procedure whereas the second is our guess at an 
additional 'fitting' systematic uncertainty from seeing the effect of changing the bin-size 
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Figure 13: Reconstructing the spectrum with a clear di-lepton edge: M q u edge after cuts selection 
B. An edge is expected at 467 GeV. Standard Model background is shown in red. 



and fit interval. For events in which M™ ax > Mu > MjJ lax /^/2, the larger qll mass is 
plotted in figure 14. 14 This gives a threshold which has a theoretical locus in terms of 



masses of all four sparticles involved in the chain [65]. The qll threshold is obtained by 
fitting with the empirical form 



[A(M - M^°) + B(M 



M^ 9h )}6{M 



M H 9h) 



(4.4) 



and a Gaussian smearing. We obtain M^ m = 188 ± 7.5 ± 6.6 GeV from the fit, which 
Fig. 14 shows, reproducing the shape of binned simulated data up to statistical variations. 

Events are then further selected such that one qll mass is less than and the other 
greater than M^ ax ~ 500 GeV. This identifies the jet involved in the decay chain qi — > 
qX2 I Pq — * Xi^ + ^ Q- By combining each of the leptons with this jet, we can plot the ql 
mass (figure 15) which constrains a function of the sparticle masses [65]. This endpoint is 
located in a similar fashion to the di-lepton endpoint, using a Gaussian smeared triangular 
fit. We obtained M™ ax = 353 ± 1.7 ± 3.6 GeV. 

We summarise the values found through the above fitting procedure in Table |2[ They 
differ from the expected values by 'experimental' systematic errors. The sources for these 
systematic errors can be associated with hadronic calorimeter calibration, jet energy leakage 
or the cone or cluster algorithms that reconstruct the jets. For our purposes we assume that 
these systematic errors can be removed as part of the experimental analysis. Following [64], 



14 The flavour subtraction of the Standard Model background in figures |l^ and [TB] and can be seen to 
slightly over subtract. This can be understood as an artifact of the triggers used, which have a small flavour 
asymmetry. 
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Figure 14: Reconstructing the spectrum with a clear di-lepton edge: heavier M q u in 4-jet events 
wither and Mu > M™ ax / 'y/2. A low-energy threshold is expected. Standard Model background is 
shown in red. 



for our estimation of sparticle masses we shift the measured endpoints to the theoretical 
values, thus eliminating this experimental systematic bias, but keep the errors arising from 
the fitting procedure. We add the statistical and 'fitting' systematic errors in quadrature 
in order to quote a total uncertainty on the central theoretical value. 

We now use these shifted end- 
points to reconstruct the XiiQL,eR 
and X2 m &sses. To do so, we take 
m q L > men and m^o to be randomly 
generated with a uniform distribu- 
tion within 50% of their central val- 
ues and compute m^o using the MJ^ ax 
di-lepton endpoint which has a very 

small statistical error. This is equivalent to approximating the very narrow Gaussian like- 
lihood distribution of M™ ax with a 5 function. We then compute the x 2 f° r the remaining 
observables (MJP q ax , M z ™ n , MJ^ ax ) and assign a weight oc e~ x ' 1 1 2 to this set of randomly 
generated masses. Doing this many times provides a sampling of a probability distribution 
for the sparticle masses. The marginalisations to three independent mass differences is 





Fitted value/GeV 


Shifted value/GeV 


M rnax 

it /[max 

1V1 llq 

Mrnax 

yrmin 
1V1 llq 


69.4± 0.15 
450±5.5±2.4 
353±1.7±3.6 
188±7.5±6.6 


69.4±0.15 
467.6±6.0 
370.8±4.0 
202.8±10.0 



Table 2: Fitted edges and their uncertainties. 



shown in figure 1€. From each probability distribution, we estimate the mass differences. 
The histograms show the binned estimated probability density function of each mass dif- 



ference, whereas the continuous lines show the best-fit Gaussian shape. In Fig. 16a, the 
plot has been separated into two regions, each of which is fitted with a separate Gaussian. 
m^o — m^o has been fitted with two half-Gaussians of different widths, "glued" at the 

A2 -M 
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Theoretical/GeV Estimated/GeV 


IR Xl 

m^o — m^o 

X.2 Xl 


37.5 28.1±1.4, 37.5±2.3 
69.7 69.7±g;| 
567 564±26 



Table 3: Sparticle mass differences. 



m^o is not very well fitted with a Gaussian, as Fig. 16c shows. We char- 

Xi 1 — 1 

acterise the distribution instead by its mean and standard deviation. The characteristic 
double-bump structure of Fig. |i~6| a displays the existence of two solutions in mass space for 
the edge variables listed in Table ^ and results in two different possible estimated values 
for the mass difference, one at each local maximum. This ambiguity was also observed for 
the case of an LHC SPSla spectrum reconstruction in [67]. We display the estimated and 
theoretical values for the three mass differences in Table [| In our case, the 'true' peak cor- 
respond to the ql-edge given by invariant mass of the quark and the far lepton in the decay 
chain, whereas the 'wrong' peak correspond to a solution for the ql-edge given by the quark 
and the near lepton. We find the 'wrong' solution by numerical scans to be rriq L = 792.6 
GeV, m^o = 303.6 GeV, mj- = 261.8 GeV and m^o = 233.7 GeV with a resulting X 2 of 
(the solution of the correct peak also has a zero value of x 2 ; by construction). The different 
heights of the bumps in Fig. [l6| must therefore be a consequence of volume effects in the 
marginalisation procedure, since the two best-fit solutions are equally likely. One requires 
additional data in order to discriminate between the two solutions experimentally. 
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Figure 16: Mass difference samplings for the reconstructed spectrum with lOOfb 1 of data: (a) 

m i R ~ m x1 ' ( b ) m x° 2 ~ m x1 ' ( c ) m M - m x? ■ 

Despite the existence of two possible good-fit regions of mass difference space, we have 
enough information to discriminate against mSUGRA models. As mentioned before, the 
gaugino mass ratios in mSUGRA are M\ : M2 : M3 ~ 1 : 2 : 6. We also know that in 
mSUGRA, rriq L is strongly correlated with m g and generically we have rriq L < m g (this is 
guaranteed within mSUGRA by the presence of light sleptons). Therefore, in mSUGRA 
models we have m„ ~ 6m^o and 

" Al 

{m g -m^o)/(m^o - m^o) = 5, 

{m-g - m^>)/(m^o - m^o) > (m 5i - m^.)/(m^o - m^o). 
However, for the large volume models we obtain 

(m q ~ L - m^o ) /(m^o - m ± o) = 8.09 ± 0.38. (4.5) 
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from Table ||. Thus the measurements made are not compatible with the mSUGRA sce- 
nario. Further measurements, e.g. of the gluino mass and the right handed squark masses, 
would provide further evidence for discrimination against mSUGRA. Using the expected 
m q~L • m g ra tios, we could also investigate the M2 : M3 « 1 : 3 prediction of the large 
volume models. For this it would be useful to directly measure the gluino mass: to this 
end the decay channel g — ► qqn — > qqxi m ay be exploited as well as the M? 2 variable [68] . 

5. Conclusions 

We have performed a detailed study of the expected superparticle spectrum and collider 
phenomenology for large-volume string models. Our main conclusions are: 

1. The large volume models give rise to a distinctive spectrum of gaugino masses, char- 
acterised by 

Mi : M 2 : M 3 = (1.5 -» 2) : 2 : 6 (5.1) 

This can be distinguished from the ratios that appear in e.g. mirage mediation or 
mSUGRA. 

The collider phenomenology depends heavily on the mass difference between Mi and 
M2 and the slepton mass spectrum. If this is large, leading to many \\ — * X?' ^ 
events, kinematical reconstruction of the spectrum is much easier. This was discussed 
in section f|. 

2. The overall spectrum tends to be more bunched than that of a corresponding mSUGRA 
model. This can be understood by the approximate unification, prior to the inclusion 
of the effects of magnetic fluxes on the brane world-volume, of scalar and gaugino 
masses at the intermediate (fundamental) scale. There is then less energy for the 
physical masses to evolve from their theoretical boundary condition and the overall 
spectrum falls within a narrower mass range. 

This effect also occurs in models of mirage mediation, where gaugino masses are 
(accidentally) unified at the intermediate scale. 

3. More concretely we find: the LSP is mostly bino. The second neutralino is mostly 
wino and is almost degenerate with charginos. Sleptons are almost degenerate, with 
stau the lightest. The gluino is the heaviest sparticle. The ratio of the gaugino-squark 
masses is larger than that predicted by mSUGRA. 

4. We have quantified the uncertainty that appears in the weak scale spectra due to 
uncertainties in the high-energy soft terms. The incorporation of such uncertainties 
is essential in trying to make predictions for LHC signatures based on high-scale 
string constructions. 

5. We have used event generators and detector simulators to study possible signatures of 
our models. We analysed the use and limitations of certain 'counting' observables to 
contrast our models with other classes of models, especially a line through mSUGRA 



- 35 - 



parameter space (the SPSla slope) that has been well studied in the literature. We 
found that some counting observables are more useful than others but in general 
they do not provide enough information to fully distinguish the models, at least in 
simple 2-dimensional projections. It may be true that a fit of the models to the full 
parameter space of counting observables is required. 

6. We studied in detail a sample model that is quite rich in \2 ~^ Xi^ T decays. 
Accurate reconstruction of some properties of the low-energy sparticle spectrum is 
possible with 100 fb _1 of integrated luminosity. Our sample study shows that the 
large volume model can be differentiated from standard mSUGRA. 

In this work we have made progress in the process of starting with a well defined class 
of microscopic models and bringing them to the point where they may be confronted with 
potential experimental measurements at the LHC. This is a positive step in the direction 
of testing classes of models derived from string theory. This is clearly a less ambitious task 
than testing string theory in its entirety, but one that may prove more fruitful. Moduli 
stabilisation with supersymmetry breaking has allowed us to find explicit expressions for 
soft breaking terms that have a well-defined microscopic origin, which is not the case for 
the well studied standard benchmark points. Issues, such as flavour universality and extra 
CP violation, that render the generic gravity mediation scenario unrealistic and have to 
be resolved by hand in most models, can now be understood in terms of the particular 
properties of string compactifications. Furthermore we have found ways to distinguish our 
models from mSUGRA and other models, using properties that may be directly measured 
within a few years of LHC running. 

There are however several open questions we need to emphasise. First of all we have 
been working on a scenario allowed by type II constructions, in which the Standard Model 
is assumed to live on a set of D-branes localised at a particular region inside the Calabi- 
Yau manifold. The problems of moduli stabilisation and supersymmetry breaking are thus 
decoupled from the details of the Standard Model construction. This approach has the 
positive feature that moduli are stabilised in a large class of models. Other issues, such as 
the number of families, proton stability and gauge unification are more model dependent. 
In this sense our results are very robust. On the other hand they lack concreteness in the 
sense that we do not have an explicit D-brane configuration with the MSSM spectrum, 
known Yukawa couplings, etc. Finding a fully realistic model in this approach is therefore 
an open question. 

Nevertheless we have identified the main sources of uncertainty in our analysis that 
parametrise our ignorance of such a realistic model. First, we included the effects of 
magnetic fluxes, usually needed to construct chiral models on D7 branes. Even though the 
dependence of Kahler potentials and gauge couplings on the magnetic fluxes is not known, 
we were able to parametrise our ignorance in terms of the random parameters e. A second 
source of uncertainty is the spectrum itself. We know that typical quasi-realistic D-brane 
models (see [69] for a recent review) usually have extra particles beyond the MSSM and 
that the hypercharge does not have a canonical normalisation. We took into account these 
effects by varying the hypercharge normalisation and finding observables, such as the ratios 
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of gaugino masses, which do not change if there are extra fields beyond the MSSM. Finally, 
we have assumed the simplest configuration of D7 branes hosting the Standard Model in 
the sense that they are all assumed to wrap the same 4-cycle. Different configurations 
would slightly change the expressions for the soft breaking terms. The expressions would 
then depend on an extra parameter A that takes different values depending upon which 
cycles the branes wrap and the manner of their intersection. For the case studied here, 
A = 1/3 [12,13]. It would be interesting to explore the implications of other configurations 
leading to different values of A. 

The model independence of our analysis makes it easy to adapt once explicit realistic 
models are constructed that may differ from the MSSM. Furthermore, potential experi- 
mental measurements at the LHC may provide guidance on what the structure of these 
realistic models should be. Even if at the end it turns out that our models will not pass 
experimental scrutiny from the LHC, the detailed analysis made, all the way from string 
theory to LHC observables, should be a useful guide for future proposals. It is encouraging 
to have this rich interplay between theory and experiment waiting for the arrival of LHC 
results. 
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A. PGS Level 2 Triggers 

Here we list the PGS Level 2 Triggers used in our event analysis for ease of reference, which 
can be obtained from pgs_olympics.f in the PGS package [9]. If any of the following apply 
to an event, the trigger is passed and the event recorded. The PGS definition of isolation 
is somewhat involved, and we refer the interested reader to the PGS manual for further 
details. Our isolation criteria for muons are based on those used in the Chameleon [70] 
package: the ratio of the energy in a 3x3 grid surrounding the muon and the pt of the 
muon is required to be < 0.1125, and the total pt of a AR = 0.4 cone region surrounding 
(but excluding) the muon is required to be < 5 GeV. 

1. Inclusive isolated lepton 1 = e,fi pr{l) > 180 GeV; 

2. For a lepton p T {l) > 130 GeV plus a jet j p T (j) > 200 GeV; 
15 Formerly PPARC (R.I.P.) 
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3. Isolated same-flavour di-leptons pr{h,2) > 60 GeV; 

4. Di-leptons pr(h,2) > 45 GeV P lus J et Pr(J) > 150 GeV ! 

5. Isolated opposite-flavour di-leptons pr(h,2) > 30 GeV; 

6. Isolated lepton pr(l) > 45 GeV plus isolated tau prij) > 60 GeV; 

7. Isolated di-tau pt{ti,2) > 60 GeV; 

8. Inclusive isolated photon pr{l) > 80 GeV; 

9. Isolated di-photon j?t(7i,2) > 40 GeV; 

10. Inclusive pr > 200 GeV; 

11. Inclusive single-jet pt > 1000 GeV; 

12. Jet pr{j) > 300 GeV plus^T > 125 GeV; 

13. Acoplanar jet p T (j) > 150 GeV and pr > 80 GeV, 1 < A(f>j P £ < 2; 

14. Acoplanar dijets pr(ji,2) > 400 GeV, A<f>jj < 2. 
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